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CN 1 Abstract 

Symmetries of nonlinear control systems in state representation are considered. To this 
end, a geometric approach to ordinary differential equations is advocated. Invariant feed- 
' back laws for systems with Lie symmetries, i.e. feedback laws that preserve the symmetry 

group of a considered plant, can be constructed based on invariants of the considered group 
action. Under minor technical assumptions suitable invariant tracking errors can be deter- 
mined by following a normalization procedure. The underlying local geometric meaning of 
this procedure is discussed and it is shown how it can also be applied in order to derive a 
local, reduced-order system representation. Further, the idea of controlled symmetries, i.e. 
imposing desired symmetry properties on a given control system by state feedback, is dis- 
cussed by outlining an exemplary control design for a predator-prey bioreactor. Symmetry, 
Lie-Backlund equivalence, Lie group, invariant feedback control 
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Symmetries of dynamical systems have been a subject of long standing interest in the treatment 
of dynamical systems. Roughly speaking, a transformation (or a family of transformations) is a 
symmetry of a dynamical system if the transformation maps solutions to solutions. The knowl- 
edge of the admitted group of symmetries can help to obtain a structural understanding of the dy- 
namics and the underlying problem that has been modeled. For instance in classical mechanics, 
the existence of sy mmetries is closely connected with the existence of con s erved quantities and 
system reduction ( Abraham and Marsden . 19871 ). In Grizzle and Marcus! ( 1985 ) the notion of 
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symmetry is adapted to nonlinear control systems in state representation (i.e. under-determined 
ordinary differential equations with input) and structural consequences of Lie symmetries in 
the light of a possible decomposition into subsystems of lower dimension plus quadrature are 
discussed. To this end, an assumed special structure of the Lie algebra (namely Abelian, non- 
Abelian with center) of the symmetry group is exploited leading to more detailed results than 
for the general Lie algebra case, which, on the other hand, is cl osely connected to the notion 



of controlle d invariant distributions and partial symmetries (cf. iNiimeiier and van der Schaftl . 



1981 Il985l ). The interest in the decomposition of a dynamical system lies in the potential sim 



plification of the control problem (i.e. dimensional reduction) as the system equations can be 
split into one part along the group orbits and a reduced order system describing the motion in 
the orbit space ("from orbit to orbit"), whe re the latter can be rega rded as a local reduced order 
model that is in some sense complete (see Zhao and Zhang . 1992). From a structural point of 
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view symmetries - as they are intrinsic (i.e. coordinate independent) pr operties of the differen- 
tial eq uations - can serve to classify dynamical systems (see for instance iLehenkvi and Rudolph 
(2004) for second order control systems). In both cases the complete knowledge of the symme- 
try algebra is assumed leading to the problem of the computation of Lie symmetries for a given 
differential equation. For control s ys tems in state representation this has been consid ered in 
Kanatnikov and Krishchenko ( 1994 ): Chetverikov. Kanatnikov. and Krishchenko ( 20021 ). Fur- 
ther. as shown in ISchlacher. Kugi. and Zehetleitnerl (|2002l ). a symmetry-based approach yields 
results on non-accessibility and non-observability of nonlinear systems described by implicit dif- 
ferential equations (this should b e related to the role of controlled invariant distributions for 
explicit systems, cf. Isidori ( 19951 )). 

In general, symmetries are not invariant under feedback, i.e. designing feedback laws using 
the usual design methods such as feedback linearization, integrator backstepping, or sliding mode 
control can lead to feedback laws that, when applied, break admitted symmetries of the control 
problem. Hence, the question arises how feedback laws can be designed that are compatible 
with a known symmetry group. 

This observation motivates the notion of invariant feedback, i.e. feedback laws which pre- 
serve the error dynamics under the action of the admitted symmetr y . To this end, an in- 
variant err or approach has been i ntrod u ced in iRouchon and Rudolph (| 19991 ) and further ex- 
amined in iRudolph and Frohlichl ()2003l ): iRudolprJ (120031). For a discussion of invarian t con- 
trol design by exact input-output linearization see iMartin. Rouchon. and Rudolph! (|2004l ) . The 
extension of the invaria n t erro r approach to invarian t asymptotic observers can be found in 
Aghannan and Rouchon ( 2002 ); Bonnabel et al. ( 20081 ). Moreover, the availability of inputs al- 



lows the introduction of symmetries into a given control problem (certain restrictions naturally 
a Pply)- Whereas this is quite naturally part of control design e.g. when compensating gravity 
effects to render a problem symmetric, the i njection of problem spec ific symmetries by feedback 
has been denoted as controlled symmetry in lSpong and Build (|2005h . In the present paper this 
approach is adopted in order to design a feedback law which is invariant w.r.t. to certain changes 
in the realization of the plant model. 

The amount of literature on symmetries of dynamical systems, especially those carrying the 
structure of Lie groups, is by far too extensive to be covered here. Details on t h e appl i cation 
of Li e groups to differen tial equations can be found in lBluman and Kumei (|l989l ): lOlverl (|l993l . 
19951); Ibragimov ( 19941 ) and the references therein. 

The intention of the present paper is twofold. On the one hand it aims to advocate a 
differential geometric approach to symm e tries based on jets and prol ongations follow ing the 
school of Vinogradov ( Vinogradov . 1981 . 19841 ) and the approach of Zharinov ( 1992 ). This 
framework allows an intuitive understanding of symmetries directly motivated from classical 
geometry, i.e. as automorphisms of some geometric object. Further, as this framework has 
been proved useful w.r. t. control-related questio ns such as the understanding of equivalence and 



differential flatness (cf. iFliess et all . 1 19941 . Il999l ). it also allows the identification of similarities 



between other geometric results. 

On the other hand, the existence of symmetries raises the questions of "suitable" tracking 
errors in order to preserve this particular structure. Whereas constructive answers to the ques- 
tion of how to design compatible feed backs w.r.t. to known Lie symmetries have been given in 



Martin. Rouchon. and Rudolph! (|2004l ). an example is given in the second half of the paper to 



point out the potential use of state-symmetries injected by suitable feedback in order to achieve 
invariance w.r.t. certain transformations. 

The article is organized as follows. In Section [2] the geometric approach to differential equa- 
tions and symmetries is outlined, followed by an intuitive definition of symmetries for differential 
equations given in Section [3l By focusing on symmetry transformations that are elements of 
a connected Lie group acting locally effectively on a manifold, invariants of these symmetries 
can be constructed by following a normalization procedure. These invariants play a prominent 
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role in the invariant tracking control approach based on so-called G-compatible tracking errors, 
which is motivated and outlined in Section 01 Further, it is shown how the normalization pro- 
cedure can be applied to derive a local reduced order realization of a control system in state 
representation admitting a state symmetry. Finally, Section [5] is devoted to a short example 
illustrating the use of symmetries induced by suitable state feedback allowing a predator-prey 
bioreactor to be rendered invariant w.r.t. different growth kinetic models. 



2 A geometrical setting for ordinary differential equations 



Since the application of differential geometric objects such as m anifolds, distributions etc. has 
become fairly stan dard in the nonlinear systems literature (see e.g. lNiimeijer and van der Schaftl . 
199(1 llsidoril . ll995l ) common objects such as manifolds and vector fields will not be define d . Gen - 
eral introductio ns to finite-dimensional differential geometry can be found in IWarnerl (11 983); 
Boothbvl <|2noah . The presentation of the geometric framework follows the setting of [Z harinov 
( 1992) which has bee n introduced in a control conte xt before, see Fliess et al. J 19941 ). Fliess et al. 
(jl999l ). IPometl (Il995l ). iBocharov et all <jl99flh . and Ida Silva. Batista, and Silverial (12007 ) for ao- 
plications, for details on the geometry of jet bundles the reader is referred to Saunders ( 19891 ). 



2.1 An introductive example 

Consider the equation of a harmonic oscillator with resonance frequency to 

Cp(t) = -uj 2 p{t), t e I C R, <p(t) 6l, (2.1) 

where the trivial bundle (M,ir = proj 1; /), M = I x M, with global coordinates (t, cp) is used 
to describe the evolution of its configuration w.r.t. time. Given a smooth function a : I —> M, 
t h-> (t,a(t)) (i.e. a graph), a is called a solution of the differential equation if a together with 
its second time derivative fulfills (|2,ip . Turning to the geometric description a solution can be 
interpreted as a graph given by a and its time derivatives in the extended space J 2 it, the second 
jet space, j 2 a : / — > J 2 ir, t h-> (t,a(t),a(t),a(t)), where j 2 a is called the second prolongation 
of a. In the jet space with local coordinates (t, cp, (p, cp) the differential equation (|2,ip defines a 
regular submanifold 

J 2 vr D S = {(t,(p,<p,<p) G J 2 vr : <p + oj 2 (p = 0} (2.2) 

introducing a geometrical object for the differential equation. As the differential equation is time- 
invariant, it is sufficient to look at the phase space for visualization purposes obtaining the well- 
known phase portrait shown in Figure [TJ A function a is a solution iff the graph of its prolonga- 
tion 

t h-> (t, a(t), &(t), a(t)) is contained in the regular submanifold: j 2 cr(t) C S, Vi € /. 

By means of the geometric approach, i.e. regarding solutions as graphs of smooth functions 
on a manifold, the understanding of symmetry of a differential equation can be related to 
the classical geometric idea of symmetry, i.e. an automorphism of a given geometric object is 
understood as symmetry hereof (e.g. the reflection about the center of a circle). Returning to 
the example, consider the mapping $ : (i, tp, cp, Cp) i-> (t, \ip, Xip, XCp), A € M\{0}. Applying $ to 
the definition of S given in (|2.2j) yields 

$(S) = {(t,p,p,Cp) e J 2 tt, : \(Cp + uj 2 ip) =0} =S, 

i.e. the submanifold is invariant w.r.t. <£. Each solution a is mapped to another solution a as S 
is invariant, and due to this characteristic property $ is denoted a symmetry of the differential 
equation. Of course this can can be observed in coordinates as well: Let Cp, Cp, Cp denote the 
transformed variables, then the transformed differential equation reads 

XCp(t) = Cp = —uj 2 \ip = —uj 2 Cp, 

showing that the differential equation is invariant w.r.t. <3?. 
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Figure 1: Differential equation (|2.1[) as submanifold S of the jet space J 2 tt (modulo time repa- 
rameterization) and action of mapping solutions to solutions 

2.2 Notions from differential geometry 

As sketched in the introductory example the geometric treatment of ordinary differential equa- 
tions in this paper aims to allow a definition of symmetry that is strongly motivated by the usual 
geometric understanding of symmetries. Hence a possible identification of a system of ordinary 
differential equations (in state space representation) with submanifolds of jet manifolds is out- 
lined based on the identification of the differential equation with its totality of solutions which 
are locally given as graphs of smooth functions. Before this geometric picture can be established 
in a more formal way some basic notions from differential geometry are recalled. 

A smooth fibered manifold is a triple (£, tt, B) consisting of the total manifold £, the base 
manifold B, and a smooth surjective submersions tt : £ — > B. The following considerations are 
restrained to the finite dimensional case, more specifically d\m.£ = 1+q, and adapted coordinates 
on £ and B are given by (t, z 1 , z 2 , . . . , z q ) = (t, z) and (t), respectively. 

loc 

For each t £ B the total manifold is locally diffeomorphic to the product space £ ~ B x Tt, 
where Tt denotes the fiber Tt = vr _1 (f) over t. Assuming that all fibers Tt are diffeomorphic to 
a typical fiber T one arrives at a fiber bundle (£, tt, B). In order to simplify notation the bundle 
is usualle only denoted by tt instead of (£, tt,B). A section a in tt is a map a : B — >• £ such that 
7r o a = id holds wherever a is defined. The set of all sections in tt is denoted by r(7r), the set 
of local sections defined around some p £ £ is denoted by r p (7r). 

In order to be able to discuss differential equations geometrically the time derivatives of 
the dependent coordinates z % have to be included within the geometric picture leading to the 
manifold of fc-jets. Two local sections a, ip E r p (7r) with cr(p) = ip(p) and further coinciding 
in their derivatives up to the fe-th order are called fc-equivalent in p. The equivalence class 
formed by all sections being ^-equivalent in p to which a belongs is called the /c-jet of a in 
p and is denoted by jlo. The reader can think of a fc-jet locally as the equivalence class of 
Taylor series expansions of smooth functions that coincide up to the A:-th order. The set of 
all fc-jets for all a € r(-7r) forms the k-th jet manifold J k ir. Adapted coordinates (t, z) induce 
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Sometimes tt is referred to as projection which it is only in adapted coordinates. 
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adapted coordinates on J k ir given by (t, z l ,Zj), i = l,...,q, j = l,...,k. Together with B 
one defines the bundle (j k 7r,irk,B). A section a € r p (7r) induces a mapping on J k ir which is 
given by its prolongation j k a : t >-> (t, cr(t), &(t), . . . , a^ k \t)). Given a local section xfj £ Tp(irk), 
ip : B — > J k ir, t h-> (t,ip(t),ifji(t), . . . ,ipk(t)) does not necessarily have to be the prolongation of 
any section in n. Since solutions of differential equations are prolongations of sections in n the 
contact distribution 

C k :u) = dz) -z i j+1 dt = 0, i = l,...,q, j = 0,1, . . . ,k - 1, (2.3) 

defined via (qk) so-called contact forms, is used to distinguish arbitrary sections in J k ir from 
prolongations of sections in it. A local section tp in J k ir is the prolongation of some section 
(j) in 7r iff the contact conditions (|2.3p hold for ip. A smooth mapping $ : J k 7r -> J k 7r that 
preserves the contact distribution, i.e. $*(C fc ) C C k , is called contact transformation, i.e. $ 
maps prolongations of smooth sections (graphs) to prolongations of (other) smooth sections, 
where denotes th e push-forward along <£. For more details on the geometry of jet spaces the 
reader is referred to ISaunderd |l98^). 



2.2.1 Differential manifold and Lie-Backlund mappings 

In the following a differential manifold is distinguished from the usual notion of smooth manifolds 
by defining a differential manifold to be a pair A4 = (M,Cm) consisting of a sm ooth manifold 



M and an involutive Cartan distribution Cm (for details cf. IZharinovl . [1992J). The Cartan 
distribution is finite dimensional and its dimension m is called the Cartan dimension of the 
differential manifold. The underlying manifold M can be in general infinite dimensional with 
coordinates indexed w.r.t. a countable index set A with cardinality ^A. If not stated otherwise 
it is assumed that each index set is finite dimensional. Note that every finite-dimensional smooth 
manifold can be made into a differential manifold in this sense by assigning the tangent bundle 
as the Cartan distribution, i.e. M. = (M, TM). A (standard) chart on A4 is given by a quadruple 
(U,ip, M m+A , Dm) consisting of a chart (U,cp) for M and a basis Dm of the Cartan distribution 
spanned by m base fields. Since the Cartan distribution is involutive for each point p € M there 
exists a unique integral manifold of dimension m passing through p. The standard charts on A4 
introduce a natural separation between the m independent coordinates (usually denoted by x l , 
i = 1, . . . ,m) and the dependent coordinates (usually denoted u a , a € A) together forming the 
adapted coordinates (x l , u a ) on M. Integral manifolds of Cm will be of particular interest as they 
will play the role of solutions of a system of differential equations. Every integral submanifold 
S C A4 can be locally written in parametric form A4 D S : u a = s a (x), a £ A, i.e. as a graph 
of smooth functions s a . In the following all considerations are restricted to the case m = 1 
(i.e. ordinary differential equations) and #A = q. The independent coordinate will be denoted 
t (time) and in order to distinguish from the general case the adapted coordinates are denoted 
by (t, z l , z 2 , . . . , z qS ) = (t, z). The set of smooth functions on M is C°°(M) and the set of vector 
fields on M with the same domain O will be denoted by T{0). 

The motivation for the distinction between smooth and differential manifolds by introducing 
the Cartan distribution can be understood from the leading example: Let us regard the plane 
S in Figure [1] as a geometric object that is to be identified with the differential equation as a 
subset of J 1 ^ with adapted coordinates (t, if, ifi) rather than in J 2 7r, i.e. elements of S that are 
prolongations of graphs that fulfill the differential equation. In general, S also contains graphs 
of prolongations that do not comply with the differential equation. The distinguishing element 
between solutions and non-solutions is their tangent space in each point. Whereas for solutions 
the tangent space is uniquely determined by the differential equation, it is not restricted for 
general graphs. If in the leading example the Cartan distribution for M = J 1 ^ is defined as 
C(M) = span{9( —co 2 ipd l p^ a differential manifold (S,C(M)) is obtained in the above sense 
consisting only of prolonged graphs fulfilling the differential equation. 
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Now, consider a smooth mapping $ : M — > N between two differential manifolds A4 = 
(M,Cm) and N = (N,Cn). The mapping is called a Lie-Backlund mapping if it preserves the 
Cartan distribution, i.e. if 3>*(Cm) C Cn holds. Regarding the coordinate form of a Lie-Backlund 
mapping the following result is recalled. 

Proposition 1 ( Zharinov . 19921 ). Let Ai = (M,Cm) and N = (N,Cn) be differential manifolds 
of Cartan dimension one. Let (U,(p,M. x M A ,5t) and (N,ip,R x M B ,U t -) on M and N denote 
their standard charts with coordinates (t,u), (t,v) and let their Cartan distributions be spanned 
by the two vector fields 



d t = d t + Y,A a (t,u)d ua , 



A smooth mapping f : U — > V given in coordinates by 



t = T(t, u) 



pel 



/3g 



(2.4) 



is a Lie-Backlund mapping iff the coordinate functions fulfill the following (defining) equations 

d t {V p ) - BP(T, V)d t (T) = 0, p G B. (2.5) 
In the following, tangency of vector fields w.r.t. submanifolds will be of importance. 

Proposition 2 (cf. Zharinovl . 1992 ). A vector field v € T(M) is tangent to a submanifold 
ScM, i.e. v\ s eT(S), iffv:T(S)^T(S), that is iff v(s)\ s = for s\ s = 0. 

A mapping / is said to be a Lie-Backlund isomorphism from A4 to M at the pair of points 
(p, q) with p S M and q € N, if / is a Lie-Backlund mapping with /*(Cm) = Cjy, and if / has a 
smooth inverse g from a neighborhood of q in iV to a neighborhood of p in M with g*(CN) = Cm- 
If this is the case, the differential manifolds M and are said to be (differentially) equivalent. 
If $ : A4 — > N is a Lie-Backlund isomorphism between two differential manifolds A4 and A^ with 
smooth inverse ^ : M — > A4, and S : N — > M is an automorphism on M, i.e. 5*(C(A/")) C C(N), 
then the mapping $oSof _1 :A4— ^Al defines an automorphism on Ai. 



2.3 Geometric interpretation of control systems in classical state represen- 
tation 

Consider the smooth manifolds £>, £ , and U referred to as base, state, and input manifold with 



adapted coordinates (t), (t, x , ; 



*) = (t,x), and (i, it 1 , . . . , u m ) = (t,u), respectively. 



From these manifolds the bundles (£, ir, B) and (U, tt u , B) are defined. Using the bundle product 
over the same base manifold x g yields the manifold M = £ XqIA being the total manifold of the 
bundle (M, ttm = n ^B^uiB) with adapted coordinates (t,x,u). The mapping p := ir* (ir u ) : 
M — > £ can be used to pull back the tangent bundle (T£,T£,£) to M yielding the pull-back 
bundle (p*(T£ ), p*(re), M) with p* (T£) ~ (t,x,u,i,x). The following commutative diagram 
summarizes the bundle construction. 



p* (T£) <- 



f 



T£ 



£ x B U < > 



£ 
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Now, consider a system of ordinary differential equations in classical state representation on 
M, i.e. 

x\t) = jx\t) = f (t, x(t),u(t)) , i = 1, 2, . . . , n. (2.6) 

The smooth functions f % 6 C°°(M) define a vector fielcd 

r(p*(T£))3v f = d t + f i (t,x,u)d x i (2.7) 

which is a section in the pull-back bundle p* (T£). A local section a : B — > M, t h4 (i, x(t),u(t)) € 
7r7j (i) whose first prolongation fulfills the system equations (|2.6p is called a solution of the 
system. Since / = (/*) depends on the input u the system is under-determined. One arrives at 
the following two possible geometric objects that can be understood as representatives of the 
differential equation ()2.6j) : 



• Differential manifold perspective: The vector field (|2.7p spans the Cartan distribution 
Cf = sp&n{vf} of the differential manifold A4/ = (M,Cf C p* (T£)), where the Cartan 
distribution distinguishes prolongations of solutions from prolongations of other sections. 
The geometric object identified with the differential equation is given by the differential 
manifold M f above. 

• Jet manifold perspective: The first jet manifold J ttm has adapted coordinates (t,x,x,u), 
and herein the differential equation (|2.6p defines a regular submanifold S C J 1 vtm given 
by 

S = {p e J l ir M | F(p) = 0} with F(p) = x - f (t, x, u) . (2.8) 

Each smooth section a : t i— > (x(t),it(i)) has a prolongation j 1 ^ to J 1 ^/ of the form 
1 1 y (t,x(t),u(t),x(t),u(t)), where a is a solution of ()2.6[) iff i x a C S holds. As geometric 
object the embedded closed submanifold S C J ttm is identified with the system (|2.6p . 

Remark 1. As mentioned before, as long as the tangent bundle is not restricted to the contact 
distribution, the submanifold S also contains sections that are not prolongations of a section in 
ttm- Consequently, it is natural to equip «S with a Cartan distribution being exactly determined 
by the contact forms ( Bocharov et al , 19991 ). thus, in the following S stands for the pair (5, TiSfl 



C 1 ). However, since in this case the Cartan distribution is not involutive, the pair does not define 
a differential manifold in the sense of the definition above, even though the construction serves 
a similar purpose. 

2.3.1 Equivalence under Lie-Backlund mappings 

By identifying differential equations with differential manifolds one obtains a notion of equiv- 
alence of systems under Lie-Backlund mappings. The following not i on of equivalence has 
been already discussed in a state space context in iFliess et all (|1994l . 1 19991 ). Two systems 



M = (M, span{i>/}) and M = (N, sp&n{v g }) are (orbitally) equivalent at (p, q) € M x N, iff 
there exists a Lie-Backlund isomorphism $ between open neighborhoods of p and q such that 
q = <3?(p). Since both differential manifolds consist of integral curves of their Cartan distribution 
the trajectories of both systems are ^-related in (p,q). The systems are (orbitally) equivalent if 
there exists such a mapping from an open dense subset U C M to N. The systems are (orbitally) 
equivalent if such a mapping exists for any pair (p, <3?(p)) on a open dense subset of M. 

Clearly, $ maps Vf to an element 4> • v gi 4> € C°°(N), of sp&n{v g }. If the parametrization of 
the integral curves is preserved, i.e. eft = 1, the two systems are called differentially equivalent. 



Throughout the paper summation over repeated indices is assumed. 
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The latter case is related to the usual notion of fli-related flows of vector fields in the finite- 



(Warner. 1983). 


Fliess et al. 


( 


1999 



isomorphisms but the state dimension n can be lost. Hence such equivalences can be useful in 
order to simplify control related problems, i.e. feedback design or the search for symmetries by 
reduction of the state dimension. 

Remark 2. It should be pointed out that it can be useful to look at a system representation of 
lower state dimension that is Lie-Backlund equivalent to the original system representation in 
order to simplify the computation of admitted symmetries. For instance, consider the equations 



y 



u 1 sin 8 + eu 2 cos ( 



u cos 9 + eu sin 9 + g, 



u 



descr ibing the planar motion of a rigid body modelling a PVTOL vehicle (cf. lHauser et al. 
1992). Here (y l ,y 2 ) is the position of the center of mass, 9 describes the orientation of the 
body, g denotes the gravitational acceleration, the inputs u 1 and u 2 are the linear and the 
rotational acceleration, and e <C 1 is a parameter describing the coupling of the equations w.r.t. 
u 2 . This system is Lie-Backlund equivalent to 



-v 2 sin?; 1 , 



v 2 cos v 1 + g, 



with z 1 = y 1 — esin#, v 1 = 9, z 2 = y 2 + ecos#, v 2 = u 1 — e9 2 . Any symmetry of the latter 
system is reflected by a symmetry of the original equations. Hence in terms of the computation 
of symmetries it might be easier to consider the representation in (z, i>)-coordinates. 

Remark 3. In Fliess et al. (| 19991 ) the notion of differential flatness of nonlinear control systems 
is presented within a differential geometric framework based on Lie-Backlund equivalence. To 
this end, consider an infinite-dimensional manifold with coordinates (?/£,), i = 1, 



, m, 



v > 0, together with the Cartan distribution C = span|^^ >0 where y l u denotes 

the is-th time derivative coordinate of the i-th component of y, i.e. a chain of m independent 
integrators of arbitrary length. A control system A4f = (M, span{i>j}) is differentially flat iff it 
is Lie-Backlund equivalent to a trivial system. 



3 Classical symmetries and symmetry groups 



As motivated in the leading example the differential geometric framework the definition of 
symmetry for differential equations can be understood from the classical geometrical meaning. 

Definition 1. A symmetry of a differential equation S C J 1 ttm is an automorphism on S. 

As pointed out in Remark [1] it is assumed that S consists of integral curves of the contact 
distribution. Therefore, any symmetry transformation is necessarily a contact transformation 
on J "Km which is equivalent to a Lie-Backlund transformation on (<S,T5 DC 1 ). 

Definition 2. A contact transformation <3? : J ttm — > J "Km for which ^(5) C S holds is called 
a classical symmetry of the differential equation S. A smooth vector field v§ G Y(TJ 1 -k) is called 
an infinitesimal symmetry of S if its flow is a classical symmetry. 

From this definition it is not clear how a contact transformation on J ttm is related to 
point transformations on M, which have been initial ly consider e d. In the finite-dimensional 
case this question is answered by a result given bv iBacklundl ()1875l ). For details see also 
Anderson and Ibragimov ( 1997 ); Bocharov et al. ( 19991 ). 
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Theorem 1 ( Backhand! . Il875l ). Any contact transformation on the jet manifold J k iTM, k > 1, 



is 



1. dim7r = q = 1: the (k — l)-th prolongation of some contact transformation on J^t^m, 

2. dim7r = q > 1: the k-th prolongation of some point transformation (diffeomorphism) of 
the independent and independent variables on J°7Tm = ^m- 

Therefore, as, n + m will always exceed one, all symmetries derive as prolongations of diffeomor- 
phisms on -km- 

Remark 4. Following the nomenclature of lBocharov et al. (jl999h the prefix classical is used for 
symmetries that preserve the order of the differential equation (i.e. the considered transforma- 
tions close off w.r.t. the order of the highest derivatives). The mentioned result of Backlund is 
exactly due to this restriction. Since this type of symmetries is defined independently from the 
differential equation, they are referred to as external - they are also defined "outside" of S. Sym- 
metry transformations that are only defined on S are consequently denoted internal, and any 
external symmetry restricted to S generates an internal symmetry. However, if one allows the 
considered transformations only to be declared on the domain of the differential equation, i.e. one 
can use the differential equation in order to re-express higher-order time derivatives occurring 
due to transformation that do not prese rve the order, one arrives at so-called dynamical sym- 
metries, which are by definition internal ( Olver. Anderson, and Kamran . 1993). More generally, 
if an infinite-dimensional framework is used (i.e. for partial differential equ ations), B acklund 1 s 
result no longer hold s allow i ng th e introduction of generalized symmetries ([Vinogradov! . 1 19841 ; 
Bocharov etal\ . Il999l : lOlveil 1 19931 ) , i.e. Lie-Backlund transformations which do not preserve the 
order of the differential equation. For ordi nary differential equations generalized symmetry can 
be related to a dynamical symmetries, see ( Olver. Anderson, and Kamran . 19931 ) for details. 



Example 1: Consider a system of ordinary differential equations in state space representation 

z l = f{t,z), i = 1,2,..., q, z e Z c R q , t e / c E, 

denned on the differential manifold A4 / with global coordinates (t, z) and the Cartan field dt = dt 
f l [t, z)d z i. An automorphism g : M / — > M. / given by 

g : (t, z) i ^ (t, z) = (9(t, z), 0, z)), (t, z), (i, z) e I x Z, 

is a Lie-Backlund mapping iff the smooth functions 8 and C fulfill the defining equations (|2.5p : 

d t C - f (t, z) d t 6 = 0, i = l,2,...,q. 

This means nothing else but 



r (t,z) 



at + gzf z ) 



dz l 
dt 



1,2,.. 



,q, 



dt dzi 

i.e. the differential equation is invariant w.r.t. any Lie-Backlund mapping on Al/. 



The definitions given above lack constructiveness in the sense that finding all symmetries 
means to solve the impossible problem of finding all Lie-Backlund mappings of the form (|2.4p that 
preserve the Cartan distribution for a given system of differential equations, or equivalently, all 
contact transformations on J 1 itm that preserve the submanifold S. However, it is well known 
that the situation eases if more structure within the family of symmetry transformations is 
assumed, namely to consider only transformations that are elements of a Lie grourjfl acting 
smoothly on S. In order to follow this approach the next section recalls some necessary facts 
about L i e gro u ps. Fo r deta ils on Lie groups and their application to differential equations see 
Warner! (| 19831 ): lOlver! (|l993l b 



Note that the additionai structure narrows the considered class of symmetries down, e.g. discrete symmetries 
like reflection are excluded. 
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3.1 Lie symmetry groups and invariants 

An r-parameter Lie group is a group G which also carries the structure of an r-dimensional 
smooth manifold in such a way that both the group operation m : G x G — >• G, m(g, h) = g ■ h, 
g,h G G, and its inverse map i : G — > G, i(g) = g^ 1 , g G G, are smooth maps. Since the group 
operation carries a smooth inverse it defines a diffeomorphism on the manifold it acts on. From 
now on, only local groups acting on a smooth fibered (1 + q)-dimensional manifold (M,n,B), 
M ~ (t, x , ... , x q ) = z, B ~ (t), are considered. 

Definition 3. Let G be a Lie group with identity e and W C M an open set. A local transfor- 
mation group {cj)g)g & G on IF is a smooth map (g, z) G G x W h-> g (z) € W such that $> e (z) = z, 
z G W, and ^g 2 (4>gi{ z )) = ^sa-siC 2 )) 9u9i £ G, z £ W where defined. 

In local coordinates one can write the action of the transformation group smoothly parametrized 
w.r.t. its r group parameters a k , k = 1,2, ... ,r, 

4> g : W x W -> W, z 0(z;a), G3j'"a= (a 1 , a 2 , . . . , a r ) G M r . 

The action of all group elements on a given point defines the orbit of a transformation group, 
i.e. a minimal nonempty invariant subset O C W . If the transformation group acts globally on 
M, the orbit O z through z G M is the set O z = {g ■ z\g G (?}. From the definition of an orbit it 
follows that a subset W C M is G-invariant iff it is the union of orbits. Clearly, points belonging 
to the same orbit form an equivalence class w.r.t. their orbit membership. A transformation 
group acts locally effectively if from g ■ z = z for all z G W follows g = e. Further, a local 
transformation group acts semi-regularly on M if its orbit dimension is constant. The group 
acts regularly if the action is semi-regular and for each point z G VF there exist arbitrarily small 
neighbourhoods U z with the pro perty that each orbit thorugh z of G intersects U z in a pathwise 
connected subset@ (j01verlll993h . The group acts free on M if for any two distinct g,h G G and 



Vz G W , g ■ z h ■ z, i.e. only the identity element e has fixed points. 

The transformation (f) g = {<p l , 4> 2 , . . . , (p 1+q ) is the flow of its infinitesimal generators 



1+9 d 
i=i 



8 zh k = l,2,...,r, (3.1) 

a=0 



which are smooth vector fields spanning the tangent space of the orbits at the identity element 
(a = 0). The key observation regarding transformation groups is that they are completely de- 
scribed by their associated Lie algebra, which in the case of effective group actions is isomorphic 
to the Lie algebra spanned by its infinitesimal generators. Consequently, instead of working 
with the local action <p g it suffices to consider the r infinitesimal generators^. 

Definition 4. A local transformation group G acting on S is a symmetry group of the differential 
equation iff g ■ S C S for all g G G. 

Now, let W C M with adapted coordinates (t, z). For each element g G G with g ■ W C W 
there exists an induced group action on the first jet bundle J 1 it corresponding to the first 
prolongation pr^g. As before in the general case Theorem [1] holds, and thus for more than one 
dependent variable any contact transformation on the jet manifold is derived as prolongation of a 
point symmetry. The symmetry condition g ■ S C S is fulfilled iff the group transformations only 
shifts solutions on S, i.e. leaves S invariant. This is the case if the infinitesimal generators of the 



4 As a necessary condition each orbit of G is a regular submanifold of M. 

5 If the group action is not effective, the r vector fields span a Lie algebra that is isomorphic to the Lie algebra 
of the effectively acting quotient group G\Gw, Gw = r\ z ^wG z — {g\g ■ z — z,Vz £ W} being the global isotropy 
subgroup of G. However, as the quotient group carries the s ame transform ation properties, one can consider the 
quotient group in case that G does not act effectively on W (|Qlverl.[T995l ). 
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transformation group are tangent to S. Combining Backlund's theorem and Proposition [2] one 
arrives at the following symmetry conditions for a system in state representation and infinitesimal 
generators defined on M: 



pr 



{1) v k (F(p)) = 0, Vp G S and k = 1, 2, 



(3.2) 



It turns out that the equations (|3.2p are n e cessar y and sufficient for 4> g being a symmetry of S. 
For a proof the reader is referred to Olver ( 19931 ). 

Recalling the idea of using invariant tracking errors for the control design, invariants of the 
considered symmetry group are of interest. 

Definition 5. An invariant is a real- valued function I : M — > K for which I(g-z) = I(z), \/g € G 
and Vz E M, holds. Further, let pr^^G be the fe-th prolongation of G acting on J k 7r, k > 1. A 
differential invariant is a real- valued function / : J k ir —¥ M for which I(pr^g -p) = I(p), Vg € G 
and p € J k ir, holds. 

From this definition it follows that every invariant function is constant along the orbits of 
G. The following theorem guarantees the existence of functionally independent invariants. 



Theorem 2 (i01verl . ll995h . Let G be a Lie group acting semi-regularly on the (q+1)- dimensional 
manifold M with s-dimensional orbits. Then, in a neighborhood W of each point zq £ W C M , 
there exist (1 + q — s) functionally independent invariant functions I\{z), . . . , Ii +q - s (z). 

The existence of invariant functions follows from the Fro benius t h eorem applied to a s- 
dimensional involutive subalg ebra of the Li e algeb ra (compare Warner ( 19831 ). Theorem 1.60). 
For a complete proof see also lOvsiannikovj (|19820 . §17.3. Further, any other invariant (local) 
function J can be expressed as function of Ii(z), . . . , Ii +q ~ s (z), i.e. J = H(I±, . . . , Ii +q - s ) for 
some suitable function H. Another important fact is that for any group G acting locally effec- 
tively on M with r > dimM there exists a number 5 > such that the action of pr^G (the 
(5-th prolongation of G) is locally free on J 5 M with orbit dimension s = r, i .e. the dimen s ion of 



the orbits on the 5-th jet manifold of M is equal to the group dimension (|Ovsiannikovl . Il98l 
324.1). 



3.2 Computation of invariants by normalization 

A complete set of (1 + q — r) functionally independent i nvariants for p v^'G can be constructed 
bllowing a n oi 
Olveri (|l999T ): 



rpr 

by following a n ormalization procedure first proposed in iKillind (|1889i ). for a modern treatment 



sec 



9pr(*)</> 9 (p;a) 



da 



a=0 



r, for p € W C J s tt, i.e. the prolonged action 



1. Find 5 > 0, such that rank 
is locally free on J s ir with orbit dimension s = r. 

2. Select r components [(f) 1 , . . . , </> r ) of pi^'4> g such that the normalization equations 

ci = cf> l (p; a), c 2 = 4> 2 {p; a) , c r = 4> r (p; a) , c. G E, (3.3) 

are regular w.r.t. the group parameter o. The local solution yields an equivariant map 
a = 7 (p) called as moving frame. 

3. Finally, a complete set of functionally invariants is obtained by using the map 7 in the 
remaining (1 + q — r) transformation equations, 



Ii = <P 3+r (p)\ a=<p y j = l,...,l + q 
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Figure 2: Local foliation of the manifold M generated by the group orbits, choice of a canonical 
element for each orbit, and meaning of 7 in terms of the group action in a local rectifying chart 



The construction of the invariant functions becomes clear in local Euclidean coordinates for 
J s tt. Since the Lie algebra (possibly of the prolonged group action) is involutive, invocation of 
the Frobenius Theorem allows one to choose local rectifying coordinates, as depicted in Figure[2j 
For simplicity, assume that the group acts locally freely for k = with orbit dimension r on 
M. In a suitable coordinate chart the first r components z 1 of z = (z^z 11 ^) are the orbit 
coordinates whereas the remaining coordinates z are constant along the orbits. By solving 
the normalization equations (|3.3p a canonical element is chosen for each orbit, z = Zq. For 

loc 

each point in p £ M the resulting map 7 : M — > G ~ W defines a group element such that p 
is mapped to the associated canonical orbit element, i.e. the first r coordinates coincide with 
Zq. Consequently, the normalization introduces coordinates along the group orbits w.r.t. the 
canonical element given by 7 - it defines a repere mobile (moving frame). As the normalization 
equations are equivariant w.r.t. the group action, their solution 7 enjoys the same property, i.e. 
"i[gz) ■ (gz) = zq = j(z) • z =>■ ^y(gz) = j(z)g . Consequently, by using 7 in the remaining 
transformation equations for z in the last step of the normalization procedure, one obtains 
the (local) canonical element of the orbit, which is an invariant along the group orbits by 
construction. 



3.3 Computation of classical symmetries of nonlinear control systems 



Together with a general ansatz for an infinitesimal symmetry of (|2.6p a set of linear first order 
partial differential equations can be derived from the symmetry conditions (|3.2p in order to 
compute the admitted symmetry group of a given system. Taking Theorem [1] into account, one 
arrives at the following ansatz for an infinitesimal generator on M 

v = £(t,x,u)dt + rf(t,x,u)d x i + ip 3 (t,x,u)d u j . (3.4) 

Its first prolongation is given by (see lOlveri (|l993h for a proof fl 

v / v / (3.5) 

with C = D t rf - x*D t t, IS = D t i(P - u j D t £, 

where Dt denotes the total derivative w.r.t. t. Note that the bracket notation [k] is used 
to denote time derivatives up to the fc-th order. Defining the vectors 77 = (77 , r/ 2 , . . . , r/ n ) , 

C = (C 1 , C 2 , • • • , C n ) T > <P = (^\ </? 2 , • • • , ip m ) T and applying (|3.5|1 to f|2.8|> in combination with 
condition (|3.2p yields 

pr (1) t> (F) = D t r] - AD t £ - f x r] - f u <p = r) t + rj x f + rj u u- f (& + £ x f + £ u ii) - f x rj - f u tp = 0. 



6 The pr^'w-notation is used instead of j 1 v since, in contrast to the prolongations before, the resulting trans- 
formation is in general not a bundle-morphism, i.e. it is not fiber- preserving, which would be the case for £ = £ (t). 
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Due to its linearity in u this system splits into two parts: 

Vu - fin = 0, 

Vt +Vxf - f (& + Cxf) - f x r] - fuV = 0. 

In general, solving this system of linear partial differential equations is as difficult as solving 
the original differential equation. Hence, finding symmetries of a given nonlinear control system 
usually involves consideration of different ansatz vector fields reflecting some kind of intuition 
e.g. in terms of the dependence on certain system variables etc. 

Remark 5. Since the vector field (|3.4p generates a point transformation onXxf/ its prolongation 
to J 1 ^ is contact by construction. Hence, all transformations generated by such a vector field 
naturally form Lie-Backlund maps of the contact distribution on J 1 it formally qualifying them 
as potential symmetries. 



4 Invariant control design for control systems with Lie symme- 
try 

Consider a control system in state representation (|2.6p and assume that it admits a local trans- 
formation group (tp g x ip g ) geG 

X x U 9 (x,u) = (ip g (x),ijjg(x,u)) , geG, 

where <p g is a local diffeomorphism and ip g is regular w.r.t. u for all x G X (coordinate change 
and regular state feedback), and the equalit}0 



^(x)f(x,u)=f(c Pg (x),i> g (x.u) l (1.1) 



holds for all g G G and (x, u) G X x U. Note that (p g necessarily depends only on x since the 
time is not transformed (compare equation (|3.6p ). Further, assume that the known symmetry 
should be preserved under a feedback law that is to be designed in order to stabilize the output 

y = h(x), h(x) = (h 1 (x),h 2 (x), h m (x)) , h G C°°(X), 

along some suitably planned smooth reference trajectorjjf] t h-> (yd(t),yd{t), ...)=: y d - 

To this end, assume that the output y = h(x) has a well-defined vector relative degree 
) ' r% = n > sucn that the decoupling matrix 



r = (r , . . . , 



is generically regular on X x U (for more details see e.g. lNiimeiier and van der SchaftJ . Il99(ll : 
Isidoril . 1 19951 ). The smooth map (:l->Zc M m , defining new coordinates z = (z 1 , . . . , z m ) , 



C: z 1 = h 1 (x), z 2 = h 2 (x), z m = h m {x), 

induces a local diffeomorphism 

z) = V f K\x), j = 0,...,r i -l, i = l,...,m, 



7 Recalling Example [T] one observes that this is equivalent to (ip g x 4> g ) g! z G defining Lie-Backlund mappings on 
X x U. 

8 For notational ease y d denotes the reference trajectory including time derivatives up to the required order. 



4 INVARIANT CONTROL DESIGN 



14 



and one obtains the following system representation in (z, u)-coordinates: 

for j = 1, . . . ,r l - 2, i = l,...,m. (4.2) 



Z) - Zj +1 , 



= F l (z,u) 



Since $ defines a diffeomorphism on X its first prolongation is a Lie-Backlund mapping w.r.t. 
the contact distribution of J 1 -k. Therefore, the local symmetry group (ip g x tp g ) g& G induces the 
local symmetry group (jp g x V's) 9gG m ( z -> ^-coordinates with Tp g = $ o ip g o 3> _1 . Defining the 
usual tracking errors e % = y % — y d w.r.t. the components of the input and the reference trajectory 
and imposing a linear time-invariant error dynamics 

e\i + cf^ei^ + ■■■ + c\& + c\e { = = E i (e\ e[, ... , e%) , i = l,...,m, (4.3) 

where e*- denotes the j-th time derivative of the i-th error component, and using the model 
equations (|4.2p one obtains a feedback law u = U (x,y d ) realizing the desired error dynamics. 
However, rewriting the error dynamics in transformed (x, £t)-coordinates only results in the 
identical tracking behavior, iff each E % , i = 1, . . . ,m, is equivariant w.r.t. the local symmetry 
group, i.e. 



f p g o E l (e\ rl] ) = E l [i T p g o e l [ri] j for all g e G, i = 1, . . . , m, 

where j rl ~p g ° E i and f l p g ° e 1 ^ denote the induced action on the error dynamics and the error, 
respectively. If these equalities hold, the stabilizing feedback law U renders the system fj4.2|) 
invariant w.r.t. the induced symmetry group. Moreover, since symmetry is independent of the 
specific system representation, the symmetry of the original state representation is preserved 
under this feedback as well. 

Equivariance is clearly achieved if the error is invariant w.r.t. the group a ction motivating an 



i nvari a nt feedbac k appr o ach based on invariant tracking errors a s proposed in lRouchon and Rudolph 



(|l99flh : iRudolnhl tooi ): iMartin. R.ouchon. and Rudolphl (j2n04T ). The given feedback design as- 



sumes an induced group action on the output y, which due to the special choice of the coordinates 
was given by Tp g . In a more general setting this motivates the following definition. 



Definition 6 (jMartin. Rouchon. and Rudolphl |2004| ). An output y = h(x,u) is called G- 



compatible if there exists an induced transformation group (p g ) g& G on Y such that h(g-(x, u)) = 
Pg ■ h for all g € G. 

The special distinction of G-compatible outputs becomes necessary as in the following the 
invariant tracking error will be defined in terms of the output y only. Consequently, the group 
action has to be restricted to Y inducing the necessity of the co-distribution generated by the 
elements of the output map h = (h 1 , . . . , /i m ) T 

Q = span {dh 1 (x),dh 2 (x), ...,dh m (x)}c T*X, 

to be invariant along the group orbits, i.e. £1 is closed w.r.t. the group action, 

we!] L Vk oj G n for k = 1, 2, . . . , r , 

where L Vk denotes the Lie derivative along the fc-th infinitesimal generator of the group. This 
observation might be exploited if G-compatible outputs are not known as shown in the following 
example. 
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Example 2: (Computation of G-compatible outputs). 

Consider the well-known equations of the kinematic car (unicycle) 



cos 9 
sin 9 



= — tan</?, 



(4.4) 



with 



being the position of the rear axle center, 9 



describing the orientation of the car, and the inputs v (velocity) and ip 
(steering angle) (Figure|3|). The model equations are invariant w.r.t. the 
special Euclidean group SE(2) consisting of rotations and translations 
in the plane with infinitesimal generators 



i>i = — z 2 d z i + z 1 d z 2 + de 



v 3 = d z 2 . 




Figure 3: Kinematic car 
! (x)}. If is invariant w.r.t. 



With x = [z , z , 6) and u = (v,p) the model equations are in 
state representation. Consider a two-dimensional smooth output y = 
(h 1 (x) , h 2 (x)^ , spanning the smooth co-distribution ft = span {dh 1 ^), dh 

the group action this is necessarily also true for its annihilator = span{«;} spanned by a smooth 
vector held w (x) = w 1 (x)d x i + w 2 (x)d x 2+w 3 (x)d x 3. The functions w l , i = 1, 2, 3, derive from h using the 
conditions ( dh 1 , «;) = 0. Computing the Lie derivatives along the infinitesimal generators and re- using 
the annihilating property one obtains a set of six linear partial differential equations for h 1 and h 2 



(dh 1 , [w,v k ]) = 0, 



1,2, 



.3. 



Using different ansatz functions for h 1 , h 2 one can for instance find that any given pair of smooth function 
h % {z x ,z 2 ), i — 1,2, forms a G-compatible output, showing that y = (z 1 , z 2 ) is a valid choice. Further, 
setting h l — h^z 1 ^) yields the solutions 

h 1 = h 1 {z 1 ,9), h 2 = const.; h 1 = h^z 1 ), h 2 = h 2 (z 1 )- 

h 1 = h 1 ^), h 2 = h 2 (9); h 2 = h 2 (z\9), h 1 = ^{-h 2 ^ 1 ,9)), 



which are not suitable for control purposes. 



Remark 6 (G-compatible flat outputs). In Martin. Murray, and Rouchon ( 19971 ) the notion of 



symmetry-preserving flat outputs is given which are G-compatible flat outputs, i.e. the considered 
symmetry group induces a local group of transformations acting on the flat output. Further, 
the flat output y and its image y under SO (2) 

y = [z 1 , z 2 + i 1 ) , y = (z 1 cos a — z 2 sin a, z 1 sin a + z 2 cos a + z 1 cos a — z 2 sin a) 

for the kinematic car are given as an example for a flat output that does not preserve the 
symmetry w.r.t. SO(2) since y is clearly not transformed on Y only. However, since symmetry is 
invariant under Lie-Backlund equivalence, any symmetry transformation acting on a particular 
flat output induces symmetry transformations on any other flat output[f|. Consequently, every flat 
output is symmetry-preserving in the suitable context. For the example above this can be seen as 
follows. As recalled in Remark [3] every flat system is Lie-Backlund equivalent to a trivial system, 
i.e. an infinite chain of integrators of width m. For the well-known flat output y = (z 1 ,^ 2 ) 
one obtains the differential manifold A4 y = (Ix IRg , spanji^}), with the Cartan vector field 
v y = dt+y\d y i +y 2 d y 2 +y\d y i +y\d y 2 + • • • . This differential manifold is Lie-Backlund equivalent 
to the differential manifold A4j = (l x I 3 x span{i?y}) defined by the model equations via 
the Cartan vector field vj = dt + v cos 9d z i + v sin 6d z 2 + % tan ipdg + vd v + (pd^ + vdy + ■ ■ ■ , which 
can be seen by expressing z , z 2 , 6 and v, ip in terms of the flat output and its time derivatives, 

z x =y\ z 2 = y 2 , # = arctan4, v = ±J{y\) 2 + (yf) 2 , y? = ±arctan[Z iM^jM 



9 In fact, any contact transformation on a flat output forms a symmetry, and hence induces a symmetry 
transformation on any other flat output of the system. 
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On the other hand, the differential manifolds Ai y is Lie-Backlund equivalent to My (which is 
defined analogously), i.e. both manifolds are transformed into each other using the mappings 



$:My^M ¥ ,t=y,y=y + </!,••• and y-.My 



My, y 1 = y 1 , 



2 —2 —1 

y =y -Vi, 



$*K) = dt + y\d t + {yl + yl)dy2 + y^i + (y\ + y\)d V 2 + ■■■ 
= d t + y\dyi + (yl -y\ + y\)dyi + y\dyi + (y| -y\ + y|)% ? + ■ ■ ■ = v v . 

Hence, the symmetry transformation y = (y 1 cos a — y 2 sin o, y 1 sin a + y 2 cos a) , a € [0,27r), 
induces the transformation 

TTJOO . TTJOO . ~1 —1 



i. 2 — i..o . y = rcosa - (y z - y{)sma, y 2 = y 1 sin a + y 2 cos a - (y? - y|) sin a, ... 

which preserves the Cartan distribution, and therefore is a symmetry w.r.t. the system represen- 
tation My, showing, as claimed before, that y is symmetry-preserving as well. However, as both 
system representations are connected via an infinite Lie-Backlund isomorphism (i.e. not a dif- 
feomorphism/ point transformation), the symmetry is not classical but a generalized symmetry 
in the aforementioned sense. 

Based on the notion of G-compatible outputs invariant tracking errors can be readily defined. 

Definition 7 ( Martin. Rouchon. and Rudolphl 2004 ). A smooth mapping {y,y d ) I (y,y d ) is 
a G-compatible tracking error if y i— > I (y,y d ) is locally invertible w.r.t. y along the reference 
trajectory, if / (yd,Vd) = holds for all y d £ Y, and if / is an invariant w.r.t. the induced group 
action (p) g eG on Y. 



Proposition 3 (jMartin. Rouchon. and Rudolphl . 120041 ) . Let G be a regular, locally effective Lie 
group and y a G-compatible output. Then there locally exists a G-compatible tracking error. 

Recalling the fact that every locally effective group can be prolonged up to some large enough 
order 5 to become locally free with orbit dimension r and invoking the normalization algorithm 
on the reference trajectory, a local moving frame g = 7(2/^) is obtained. Using 7 in the remaining 
transformation equations for y yields a set of m functionally independent invariants Ij (y,y^) = 
Pg • y J \ g= ^y d y Further, invariant tracking errors are given by e J = Ij (y,y d ) - Ij (yd,y d ), 3 = 
1, . . . ,m. 

Once a complete set of invariant tracking errors has been derived, invariant feedback laws, 
i.e. feedbacks that preserve the symmetry, can be designed using well-known design approaches 
such as input-output linearization, integrator backstepping, or sliding mode control. 

Example 3: (Invariant tracking errors obtained by normalization). The following tracking errors can 
be derived by normalization: 

• kinematic car, planar motion of a rigid body: The symmetry w.r.t. SE(2) is preserved when a mov- 
ing frame attached to the reference trajectory for a flat output y is used, e.g. 

e = ({y — yd, t) , (y — yd, v)), with t H> yd(t) denoting a smooth reference trajectory with well- 

define d tangent and normal vector r a nd v, see for instance Woernlej ( 19981 ): Martin. Rouchon. and Rudolphl 
(l2004h : lR,udolDh and Frohlichl (l2003h . 



attitude control of a rigid body (e.g. satellite): Let R, Rd G SO(3) denote two rotational matrices 
describing the current and desired orientation of a rigid body w.r.t. some inertial frame. Using 
the normalization equation R g Rd = R e = I, with R g being an element of SO (3) representing the 
group action and / denoting the identity matrix, one obtains the map R g = i? 7 (# d ) = Rj yielding 
the invariant error RjR. Using a quaternion approach to para metrize the rotational matrices on e 
arrives at the well-known invariant rotational error e = q d ~ 1 q ( Wen and Kreutz-Deleadol Il991 ). 
More generally, if one deals with a system on Lie groups this construction holds generally, i.e. the 
invariant error is given by o j ■ g derived from q ■ Qd = e with e bein g the identity element of the 
Lie group (see for instance Maithripala. Berg, and Davawansa , 20061) . 



scaling group: Using a relative error when a symmetry acts as a scaling on the output yields an 



invariant error e 



1, see for instance iRouchon and Rudolphl (|199 
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4.1 Structural consequences: system representation of reduced order 

As pointed out in IGrizzle and Marcus (|l98l Il985h . the existence of Lie symmetries imposes 
consequences on the underlying structure of the cont rol system, which is closely related to the 
notion of a c ontrolled invariant distribution (compare iNijmeiier and van der Schaftl . Il98l Il985l : 
Isidoril . ll995h . In order to give an example of further application of the normalization procedure 



it is necessary to recall some facts regarding the structural consequences of Lie symmetries. 
Consider a so-called state- symmetry, i.e. a symmetry group acting locally freely only on the 
state variables with orbit dimension r. Since the infinitesimal generators form a Lie algebra 
their Lie brackets w.r.t. each other can be rewritten in terms of infinitesimal generators and the 
structure constants £ 1 of the Lie group, 



cjjVi + cfjV 2 + 



+ c\jV r . 



Note that the Lie algebra forms an involutive distribution on X and hence, there exist a local 
chart for X such that the orbits of G form submanifolds with x 1 = const., i = r + 1, . . . , n, i.e. 
the group acts only along the first r coordinate directions. The construction of such a chart is 
for instance given in IWarnerl (|l983h along with the proof of Theorem 1.60. In this new chart 
the infinitesimal generators have the form 



Vi = ip} (x 1 ,.. .,x r ) d x i + IP- (x 1 , . . .,x r ) d x 2 H h (fifa 1 , ...,x r ) cV, i 

Consequently, the symmetry conditions pr^Vi(x — f(x,u)) on S read 



1 r 

5***5 





M 


dx 1 


dx r 


M 




dx 1 


dx r 



V 







(n-r)xr 







rx(n-r) 







(n—r)x(n—r) ) 



(f X \ 

f 



dx 







[vf,Vi] = 0, 



(4.5) 



for i = 1, . . . , r, i.e. the infinitesimal generators commute with the vector field Vf = f l d x i. From 
this it follows that the Lie algebra q = span{t;i, . . . v r } is invariant along the flow of Vf, i.e. 

€ for all v € - it forms a controlled invariant distribution of the control system. 
Examining the Lie brackets of the elements 3 ej = d x i, i = 1, . . . , r, with the vector fields Vf 
yields the conditions 



e; 



dx 1 



dx 1 



dfl 

dx 1 



Xj(x\...,x r )d x i + --- + X r i 



")d x r 



which have to be fulfilled for suitable A^. Hence, the functions / r+1 ,...,/ n have the form 



.r+l 



, . . . ,X",U 



arrives at the system representation 



, u m ) , k = r + 1, . . . j n in this particular coordinate chart, i.e. one 
with£= (x r+ \...,x n ) , 77= (x 1 , 



,x 



From this local representation a possible decomposition of the control system into two parts can 
be observed: 

• a (n — r)-dimensional system £ = F\(£,u) describing the motion transverse to the orbits 
(i.e. "from orbit to orbit"), and 



an r-dimensional system f\ = F2(£,r],u) describing the motion along the orbiti 



10 



10 Since r\ locally defines coordinates on G, i.e. g ~ rj, this part of the motion can be interpreted as motion on 
the Lie group. In fact, if the decomposition exists globally one can express this part as dynamical system on G 
as mentioned in Remark [7] 
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Depending on the control problem, it may be sufficient to consider the (n — r)-dimensional 
reduced order system, as the motion along the orbits is completely det ermined by the evolutio n 
of £ (and possibly some remaining components of the input), see also Zhao and Zhang ( 19921 ). 
The coordinate transformation yielding the (n — r) coordinates £ , . . . ,£ n-r can be constructed 
using the normalization procedure. Solving r normalization equations (possibly renumbering 
the components of the state) 



Ifg ■ X = C , <fg 

yields a moving frame g = 7(x 1 , . 
mation is obtained from setting 



c l 6 



r ) . By the usual procedure the local coordinate transfer- 



r+l I 

■ x I^ObI, 



,x r ) 



e 



=7(x 1 ,...,x r ) 



allowing to introduce the reduced order representation 

e 



(4.6) 



=(ci,..., C -,H-i(0) 



with S = (H 1 , . . . ,H n r ). An application of this procedure is given in Example [H 
Remark 7. For a global reduction it is necessary to construct a global decomposition of vj along 
and transversal to the group orbits. A suitable geometric formulation can be given by using 
the bundle (X x U,ttg,X/G x U), ttg ■ X x U — > X/G x U, itg(x,u) = (O x ,u). Assuming 
that the group acts (globally) freely and properly on X with orbit dimension r, the orbit space 
X/ G is in fact a smooth manifold and ttg is a submersion, justifying the bundle construction 
(cf. Abraham and Marsdenl ( 19871 ) . Proposition 4.1.23). Therefore, a unique decomposition of 
Vf into a horizontal vector field being a section in T{X/G) and a vertical vector field is possible 
based on the global foliation of X generated by the group orbits. Hence, the quotient system on 
X/ GxU can be defined globally. Further, the remaining r-dimensio nal systems can be rewritten 
as s ystem on the Lie group G. For details th e reader is referred to Grizzle and Marcus ( 19851 ) 
and lNiimeiier and van der Schaftl (jl98l Il98fih . 



Example 4: (Kinematic car: reduced order state representation) 
Based on the rotational invariancc of the model equations (|4.4[) 
consider the normalization equation 







7(0) 



Solving for the group parameter a (angle of rotation) and using 
the resulting map 7 in the transformation equations for y yields 
the new (rotational invariant) coordinates 



z = R 



(-e)V 



Rjy, 



where Rg G SO (2) denotes the rotational matrix w.r.t 6. 
tion (|4.6p yields the equations of motion in z-coordinates: 

cos6> sin#\ ( cos# N 
- sin 9 cos 8 J \ sin 8 



Equa- 




Figure 4: Interpretation of the z- 
coordinate frame 



smf 
cos ( 



COS f 

— sin( 



/ COS ( 

I sin£ 



■sine 



6=0 



v 

j tan ip 



This can be interpreted as follows: The new z describe the position of the rear axle midpoint w.r.t. a 
coordinate frame that is rotated about the origin by 8 and is rotating with angular velocity —8 (8 = 0). 
Since the translation along the tangent direction with velocity v and the rotation of the coordinate frame 
are decoupled motions, the motion of z is described by the superposition of both (Figure EJ. The nor- 
malization equation also delivers the equation for the motion along the group orbit. From + a = 8 one 
observes that is is described by the third model equation, i.e. g ~ a, d = j tamp. 
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Figure 5: Growth rate as function of the concentration b for the Michaelis-Menten kinetic (left) 
and for the Haldane kinetic (right) 



Application of invariant feedback design to a predator-prey 
bioreactor 



Turning to the application of invariant tracking errors, consider the model of a continuously 
operated (chemostat) bioreactor containing two populations P and B of microorganisms of con- 
centration p and b respectively. The predator population P feeds upon the prey population B 
whereas the population B lives on some substrate S of concentration s. The reactor is fed by a 
medium containing the substrate of concentration sf at the dilution rate D = y, where Q and 
V denote the flow rate and the volume of the reactor content. A simplified model describing the 
evolution of the two populations reads ( Pavlou . 19851 ) 



p = —Dp + v{b)p 

b = —Db + n(s)b — av(b)p 

s = D (sf — s) — [3[j,(s)b 



(5.1) 



where a and j3 are the yield coefficients for production of both populations, and v and \x are 
bounded, monotonically increasing functions modeling the specific growth rates of P and B. In 
the following two types of growth models are considered (Figure [5]): 

1. Incorporating asymptotic saturat ion effects occurring a t high er concentrations the well- 
known Michaelis-Menten kinetic ( Michaelis and Menten . 19131 ) is applied, given by 



z^m(s) = v n 



L + b' 



with v m > 0, L > denoting the maximal growth rate and saturation constant, respec- 
tively. 



2. If in addition to saturation also inhibition effects occur, the Haldane kinetic (jHaldand . 
193(3 ) 



Mb) 



v m b 



b + K s + Kib 2 



is used with maximal growth rate v = z^h(&) = v n 



l + 2y/KJKJ " V K i 

using the saturation constant Kg and the inhibition coefficient Kj. 



at b = \ -jfi- parametrized 



The system (|5.ip possesses three kinds of equilibrium points: washout of both populations for 
p = b = 0, s = sf, washout of the predator population for p = 0, b > 0, sf > s > 0, and 
coexistence of the two populations for p > 0, b > 0, sf > s > 0. The reactor is operated at the 
third kind of equilibria which all ows the t wo p opulations to coexist. For a detailed analysis of 
the equilibrium points, please see iPavloul ()1985l ). 
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In the following, it is assumed that the control objective is to stabilize the reactor around 
equilibrium points in order to reject unmodeled disturbances using the output y = (p, 6) and 
the input u = (D,sp). Further, the growth rate for b is modeled using the Michaelis-Menten 
kinetic, i.e. 

s 

M s ) = fa y , , K, fj, m > 0. 
K + s 

The stabilizing control should be realized independently of the growth kinetic v used to model 
the growth of P. Consequently, the feedback design intends to render the feedback invariant 
w.r.t. the change from one kinetic to the other, i.e. smooth changes of the inhibition coefficient. 
To this end, the system equations (|5.ip are augmented by 



A'/ 0. „(&) l,J ' 



b + K s + Kib 2 ' 

Changing from one growth kinetic to the other can then be interpreted as point transformation 



u m b 



b + K s + \Kj + a)b z 

acting on the inhibition coefficient Kj, where a = realizes the Haldane kinetic and a = —Kj 
realizes the Michaelis-Menten kinetic. Using the invariance condition (|4.ip and choosing p and b 
to be invariant under the model change, i.e. p = p and 6 = 6, one obtains determining equations 
for the induced symmetry transformation on the substrate concentration s, the dilution rate D, 
and the feed concentration sp: 

p (i/ H (6) ~D)=p (v(b; a)-D\ ^D = S(b, D; a), 



6 (fi(s) — D) — Q!I/h(6)p = 6 (^fi(s) — D^j — aD(b; a)p s = a(p, 6, s; a), 

^p Mb) - D) +^ (6 Gu(s) — D) — au n {b)p) + 
— (D (sp — s) — (3fj,(s)b) = D (§f — s) — /3fi(s)b =4> sp = crp(p, 6, s, sf; a). 

OS 

Consequently, in order to render the mapping 

Kj = Kj + a, P = P, b = 6, s = a(p, b, s; a), 

a symmetry of the reactor equations an induced feedback (D, sf) = (S(b, D; a), &f(p, 6, s, sf', a)) 
has to be applied^!. However, the system equations are locally feedback equivalent to a system 
for which the demanded symmetry acts only on the state variables x = (p,b,s) T . Again, this 
representation can be obtained by using the normalization algorithm. Normalizing the induced 
transformation s = a(p,b,s;a) = const, yields a moving frame a = "f(p,b,s), from which one 
derives the regular feedback (diffeomorphism w.r.t. u), 

v = V(x,u) = (5(b,D;a),a F (p,b,s,s F ;a))\ a=l{pAs) , 

with the new input v being invariant along the group orbits. As the symmetry has been chosen 
to leave the output y = (p, 6) invariant, the usual set-point error e = y — yo w.r.t. an equilibrium 
point yo = const, is a suitable choice for an invariant feedback design based on input-output 
linearization. 

Remark 8. Following the spirit of the given example it is possible to define errors that are 
invariant w.r.t. set-point changes, i.e. stabilizing feedback laws that realize an identical error 
dynamics around a predefined set of se t -point s. For a set-point invariant integrator backstepping 
design please see ICollon and Eudolphl (|2O10h . 



11 In ISpong and Build (2005) symmetries that are obtained by suitable feedback are denoted as controlled 
symmetries. 
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6 Conclusion 

The advocated geometric approach to ordinary differential equations allows an intuitive inter- 
pretation of symmetries and the notion of Lie-Backlund equivalence. As an important special 
subclass of classical symmetries Lie symmetries of nonlinear control systems have been consid- 
ered, including their computation using their infinitesimal generators. Carrying out feedback 
design based on invariant tracking errors leads to control laws preserving the symmetry group 
admitted by the consi d ered sy stem. Following this invariant t rackin g approach proposed in 
Rouchon and Rudolph ( 19991 ): Martin. Rouchon. and Rudolph the present paper con- 



tributes to the application of the normalization procedure used to construct invariant tracking 
errors by giving a geometric interpretation of the so-called moving frame in terms of the group 
orbits and the meaning of G-compatible outputs. It turns out that the normalization procedure 
can also be applied in order to determine suitable coordinate transformations for a local reduced 
order realization of a control system with Lie symmetries. Finally, the idea of controlled sym- 
metries, i.e. "injecting" symmetries into a given control problem by suitable state feedback, has 
been illustrated by an example of an invariant control of a bioreactor. 
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